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ABSTRACT 

We report CO detections in 17 out of 19 infrared ultraluminous QSO (IR 
QSO) hosts observed with the IRAM 30m telescope. The cold molecular gas 
reservoir in these objects is in a range of 0.2-2.1 xlO^^M© (adopting a CO-to-H2 
conversion factor aco = O.SMq (K kms~^ pc^)~^). We find that the molecu- 
lar gas properties of IR QSOs, such as the molecular gas mass, star formation 
efficiency (Lfir/Ivqq) and the CO (1-0) line widths, are indistinguishable from 
j those of local ultraluminous infrared galaxies (ULIRGs). A comparison of low- 

and high-redshift CO detected QSOs reveals a tight correlation between Lfir and 
-^co(i-o) QSOs. This suggests that, similar to ULIRGs, the far-infrared 

emissions of all QSOs are mainly from dust heated by star formation rather than 
by active galactic nuclei (AGNs), confirming similar findings from mid- infrared 
spectroscopic observations by Spitzer. A correlation between the AGN-associated 
bolometric luminosities and the CO line luminosities suggests that star formation 
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and AGNs draw from the same reservoir of gas and there is a hnk between star 
formation on ~ kpc scale and the central black hole accretion process on much 
smaller scales. 

Subject headings: galaxies: active — galaxies: evolution — galaxies: high- 
redshift — galaxies: starburst — radio hues: galaxies 

1. INTRODUCTION 

It has become increasingly clear that the growth of central supermassive black holes 
(SMBHs) and their host spheroids are closely related since the mass of the SMBH correlates 
well with the properties of the hot stellar component of the host galaxy, e.g. the galactic 
velocity dispersion (Mbh-o"*), the luminosity /mass of the spheroidal component of the host 
galaxy (e.g. Magorrian et al. 1998; Ferrarese & Merritt 2000; Tremaine et al. 2002; Haring 
& Rix 2004). Hydrodynamical simulations on gas- rich galaxy mergers, incorporating star 
formation, SMBH growth and feedback from both supernovae and SMBH, provide a plau- 
sible picture of how galaxy formation and the growth of SMBHs are related to each other 
(e.g. Hopkins et al. 2006). However, much of the gas physics in numerical simulations, in 
particular how the star formation and central active galactic nucleus (AGN) activities feed 
energy and momentum back into the gas, is uncertain. Such feedbacks may heat up the cold 
gas and drive a powerful galactic wind or outflow, which may limit or even terminate further 
star formation and black hole growth in galaxies. 

Hao et al. (2005, hereafter H05) studied a sample of type 1 active galactic nuclei (AGNs) 
selected from local ultraluminous infrared galaxies (ULIRGs, the 8 — 1000/xm integrated 
infrared luminosities are larger than IO^^Lq) and referred these as IR QSOs. By comparing 
IR QSOs with optically selected PG QSOs and narrow-hne Seyfert Is (NLSls) in the relation 
of the IR luminosities versus the optical luminosities, they find that IR QSOs show mid-IR 
and far-IR (especially far-IR, hereafter FIR) excess. The FIR luminosities for most IR 
QSOs are larger than half of their central AGN bolometric luminosities. In contrast, the 
median ratio of the FIR to central AGN bolometric luminosities for local classical QSOs is 
four times smaller than that for IR QSOs. Comparisons of the FIR spectral index of IR 
QSOs with those of PG QSOs indicate that the FIR excess of IR QSOs relative to the PG 
QSOs is suggestive of massive starbursts. These IR QSOs not only have massive starbursts 
occurring in their host galaxies, but their optical spectroscopic and X-ray properties also 
exhibit characteristics of young, growing QSOs with high accretion rates to their central 
black holes (Zheng et al. 2002; H05). Moreover, images of these IR QSOs show that most are 
undergoing clear tidal interactions or major mergers (Canahzo & Stockton 2001; see also Fig. 
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1). The optical spectra show obvious blueshifts in permitted emission hnes in some IR QSOs, 
indicating outflows. More recent IRAM Plateau de Bure Interferometer (PdBI) and Herschel 
observations of Mrk 231, the closest IR QSO, reveal clear outflows with velocity larger than 
1000 kms~^, as inferred from broad wings of the CO line or OH P-Cygni profile (Feruglio 
et al. 2010; Fischer et al. 2010) and the coexistence of active SMBH accretion and strong 
starburst (van der Werf et al. 2010). The mid-IR spectra of IR QSOs from Spitzer/lnhaied 
Spectrograph (IRS) also show that the typical slope of spectral energy distribution (SED) 
for IR QSOs is between those of classical PG QSOs and ULIRGs (Cao et al. 2008). All the 
properties of IR QSOs are consistent with their being at an important transition stage from 
mergers with massive starbursts (ULIRGs) to dust-enshrouded QSOs and later luminous 
QSOs, and eventually early-type galaxies (Sanders et al. 1988a; Sanders & Mirabel 1996; 
Tacconi ct al.2002). Studies of these objects may thus provide significant insights into the 
understanding of how the growths of spheroids and SMBHs are inter- connected and how the 
Mbu-M^ relation evolves as a function of rcdshift. 

The conclusion that the FIR excess in IR QSOs is primarily due to starbursts requires 
a large reservoir of cold molecular gas in IR QSO hosts. In fact, up to now only 6 local IR 
QSOs (I Zw I , Mrk 231, Mrk 1014, 3C 48, IRAS 1613+658 and PG 1700+518) have been 
observed (Solomon et al. 1997; Schinnerer, Eckart & Tacconi 1998; Krips et al. 2005; Evans 
et al. 2001, 2009). The molecular gas masses of these IR QSO hosts are few times 10^ to 
1O^°M0, assuming the same conversion factor from the CO luminosity to H2 mass as that 
for ULIRGs, i.e. aco — 0.8Mq (Kkms~^ pc^)~^ (Downes & Solomon 1998). In contrast, the 
molecular gas masses of CO detected classical QSO hosts are at most a few times 10^ Mq 
(Evans et al. 2001; Scoville et al. 2003). The larger molecular gas reservoirs in these IR QSO 
hosts mean that there are rich of fuel to sustain the star formation and may support the 
speculation by Hao ct al. (H05) that the FIR excess in IR QSOs is mainly due to starbursts. 
However. CO observations for a larger sample of IR QSOs are needed to statistically confirm 
this scenario. 

As we will show later, in the IR QSO phase the SMBH and spheroid are both growing 
rapidly. Two questions naturally arise: What is the timescale for the consumption of the 
molecular gas to sustain the massive starbursts? Can the local M^n-M^ relation be estab- 
hshed/maintained in this phase? Such a spheroid and central SMBH co-evolution process is 
expected to be more important at high redshift according to the downsizing galaxy formation 
scenario (Cowie et al. 1996). While a few tens of high- redshift QSOs have been detected in 
CO (Solomon & Vanden Bout 2005; Coppin et al. 2008; Wang et al. 2010; Riechers 2011a; 
PoUetta et al. 2011), surprisingly only half a dozen local IR QSOs have been observed before 
this work. Clearly a larger local IR QSO sample detected in CO is desirable for compar- 
isons with high-redshift CO detected QSOs. For all the above reasons, we have initiated a 
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program to observe a larger sample of IR QSOs with the IRAM 30m telescope in order to 
systematically investigate the molecular gas properties of these objects. 

This paper presents the results of our CO observations of 19 IR QSOs with the IRAM 
30m. The paper is structured as follows. In Sections 2 and 3, we describe our sample and 
outline the observations and data reduction respectively. We present our main results and 
discussion in Section 4, and a summary in Section 5. Throughout this paper we adopt a 
cosmology with a matter density parameter = 0.3, a cosmological constant = 0.7 and 
a Hubble constant of Hq = 70kms~^Mpc~^. 

2. SAMPLE 

All of our 19 newly CO observed IR QSOs are selected from the IR QSO sample of H05 
and have redshift z < 0.2 except for F13218 + 0552 (z=0.205), F18216 + 6419 (z=0.297) and 
F02054 + 0835 (z=0.345). This IR QSO sample (H05) was compiled from the ULIRGs in the 
QDOT redshift survey (Lawrence et al. 1999), the IJy ULIRG survey (Kim et al. 1998), 
and the cross-correlation of the IRAS Point Source Catalog (PSC) with the ROSAT All-sky 
Survey Catalog. We also include six previously published CO detections of IR QSOs. In 
total, our sample consists of 25 IR QSOs. As the fraction of IR QSOs in ULIRGs is less than 
10% and the PSC IRAS redshift survey (Saunders et al. 2000) galaxy sample includes about 
1000 ULIRGs, so the total number of IR QSOs expected in the northern sky is < 100. Our 
IR QSO sample includes about a quarter of the total, and thus should be a representative 
IR QSO sample in the local universe (see Zheng et al. 2002 for more detailed descriptions). 
Fig. [T] shows true color images of 12 IR QSOs (out of 25) cross-identified from DR7 of the 
Sloan Digital Sky Survey (SDSS). Most of these show tidal features and may be in the final 
merging stage. 

For comparisons, we also compiled several samples with CO detections, including lo- 
cal ULIRGs, low-redshift classical QSOs and high-redshift (sub) millimeter (mm) luminous 
QSOs. Notice that we correct all the published results to our adopted cosmology. In the 
following, we briefly summarize these samples. 

(1) A ULIRG sample taken from Solomon et al. (1997) consists of 36 objects in the 
redshift range of z = 0.03 to 0.27. The CO observations were made with the IRAM 
30m telescope. There are more ULIRGs with CO observations (e.g., Chung et al. 2009; 
Mirabel et al. 1990), but here we restrain ourselves to systems with IRAM 30m data. 

(2) A local classical QSO sample with CO detections is compiled from 11 PG QSOs (Evans 
et al. 2001; Scoville et al. 2003) and 16 HE QSOs (Bertram et al. 2007), after 
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excluding several IR QSOs. Evans et al. (2001) selected their sample by IR excess 
(L(8-1000/xm)/L(0.1-1.0/im) > 0.36) in the redshift range of 0.04 < z < 0.17, while 
Scoville et al. (2003) selected their PG QSO sample with criteria oi z < 0.1 and 
Mb < —23 mag. The HE QSO sample is a nearby low-luminosity QSO sample with 
redshift z < 0.06, drawn from the Hamburg-ESO survey for bright UV-excess QSOs 
and observed also by the IRAM 30m telescope (Bertram et al. 2007). 

(3) The (sub)mm luminous high-redshift QSO sample of 29 objects is taken from Wang et 
al. (2010). They observed six (sub)mm luminous QSOs with z~6 by the IRAM PdBI; 
for each object, either the CO(6-5) or (5-4) line emission has been detected. They also 
included two QSOs at z~6 previously detected in CO and another 21 high-redshift CO 
detected QSOs in the range of 1.4 < z < 5 from the literature. Notice that this CO 
detected, high-redshift (sub)mm luminous QSO sample includes nearly all such objects 
up to 2010. 

3. OBSERVATIONS, DATA REDUCTION and PARAMETER 

ESTIMATION 

The CO observations were carried out with the IRAM 30ni telescope on Pico Veleta 
near Granada, Spain. Although the first of such attempts was conducted remotely from 
the IRAM headquarter in Grenoble in July 2007, the weather was mostly bad then and 
essentially httle useful data were obtained. The data used here were mainly collected during 
the summer and autumn of 2008. Both the 3mm and 1mm dual-polarization receivers were 
used together with the 512x1 MHz and 256x4 MHz filter banks to simultaneously obtain 
both CO (1—0) and CO (2—1) emission lines. The weather was mostly excellent during these 
observing runs, but some CO (2—1) lines were still not observed for a few targets since the 
line frequencies were redshifted out of the tuning range of 1mm receivers. All observations 
were done with a wobbler switching mode with a throw rate of 0.5 Hz and a wobbler throw of 
60" in azimuth. System temperatures were on average Tgyg ~ 150 K and 450 K at ~ 105 GHz 
and ~ 210 GHz, respectively. The total usable on-source integration time is ~50 hours from 
the 2008 runs. Table 1 summarizes the details of CO observations for all 19 sources. 

We reduced the data using the IRAM software CLASS/GILDAS. First, we discarded 
a few bad scans with strongly distorted baselines, fixed some bad channels before summing 
up all scans to look for the hkely detected CO hne emission windows. We then subtracted 
linear baselines in each of the remaining scans with the adopted velocity windows before 
summing up to obtain the final spectra that arc ftirther subtracted with linear baselines. In 
order to increase the signal-to- noise ratio, all the spectra were Hanning smoothed to velocity 
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widths of > 20 km s~^. The amplitudes of each spectrum were converted to main beam 
temperatures, T-^h, by multiplying by the ratio between the forward efficiency and main 
beam efficiency Fe^/Beg. The main beam efficiency -Befr = 0.75 and the forward efficiency 
Feff = 0.95 were applied to at ~ 105 GHz and B^s = 0.57 and F^g = 0.91 to at ~ 
210 GHz, respectively. For point sources, flux densities were obtained using the conversion 
factor = 4.95 Jy/K for the IRAM 30m telescope. 

The CO line luminosities L'qq for IR QSOs were calculated following Solomon, Downes, 
& Radford (1992) 

L'co = 23.5 Qs^B Dl Ico (1 + z)'^ [K km s"^ pc'] , (1) 

where Iqo is the velocity-integrated line intensity in units of K km s~^, Z^l is the luminosity 
distance measured in Mpc, and fls*B is the solid angle of the source convolved with the 
telescope beam in square arc-seconds, i.e. ^^s*b=4^(^5 + ^b); where 9b = 24" at 105 Ghz. 

The CO luminosities for PC and HE QSOs were estimated using the same method and 
were provided in the respective papers described in §2. We simply adapted them to our 
adopted cosmology. For the high-redshift QSO sample collected by Wang et al. (2010), 
a different approach was applied since most of these objects were detected at high orders 
(J > 5) of CO transitions. Wang et al. (2010) calculated the -^co(i-o) values for all high-z 
sample QSOs under the non-thermalized excitation assumption and adopted the line ratios 
of CO (6-5)/CO (1-0) and CO (5-4)/CO (1-0) as 0.78 and 0.88, respectively. Since only 
the CO ground state transition CO (1—0) traces molecular gas at low excitation and provides 
a direct estimate of the total cold molecular gas mass (Riechers et al. 2007, 2011b; Gao & 
Solomon 2004a), wc will use either the observed and/or converted CO (1—0) properties for 
all QSOs at different rcdshifts in this work. 

Apart from CO line luminosities, FIR luminosities arc needed as well in this work. The 
FIR luminosities (Lpm, the luminosity in the range of 42.5-122.5//m) were calculated based 
on the flux densities at 60yum and lOOyum from the IRAS Faint Source Catalog for local 
IR QSOs, ULIRGs and PG+HE QSOs, following Helou et al. (1988). For high-redshift 
(sub)mm loud QSOs, their FIR luminosities were provided by Wang et al. (2010; private 
communication); they derived Lpm based on (sub)mm continuum measurements under the 
assumption that the rest-frame FIR SED can be described by a grey-body spectrum with 
a dust temperature of 47 K and emissivity index of 1.6 (Beelen et al. 2006) except for 
J1148+5251 and J0927+2001. For these two objects, Wang et al. (2010) used a fitted 
temperature of 56 K and 46 K respectively. 

We also need the AGN-associated bolometric luminosities Lagn for our QSO samples 
(hereafter we use -Lagn to denote the bolometric luminosities from AGN for all QSOs). 
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For local PG QSOs and IR QSOs, they were estimated from the extinction corrected con- 
tinuum emission at 5100A and adopting a bolometric correction factor of 9, i.e., Lagn ~ 
9ALa(5100)A (H05). For HE QSOs, the central AGN bolometric luminosities Lagn were 
derived from their Bj magnitudes with a bolometric correction factor of 9.74 (Vestergaard 
et al. 2004). For the 11 high-redshift QSOs in Hao et al. (2008), Lagn were taken from that 
paper directly, calculated in the same way as for HE QSOs. For the eight z ~ 6 QSOs, the 
absolute AB magnitudes at rest frame 1450A derived by Wang et al. (2010) were converted 
into the B band magnitudes according to Schmidt et al. (1995; see also Fan et al. 2001), and 
then a bolometric correction factor of 9.74 was employed to convert the B-band magnitudes 
into -Lagn- For the ten high-redshift QSOs neither in Hao et al. (2008) nor in Wang et al. 
(2010), we computed their Lagn in the same way as Hao et al. (2008). 



4. RESULTS AND DISCUSSIONS 

4.1. CO Emissions from IR QSO hosts 

17 out of the 19 new IR QSO hosts have been detected in CO (1—0), while nine of 
them have also been observed and all detected in CO (2—1). All the detections except 
for F12265+0219 and F12134-f5459 (which were marginally detected at 3a) have signal-to- 
noise ratio (S/N) greater than 5a. Fig. [2] presents the CO (1—0) and CO (2—1) spectra for 
all IR QSOs. Those targets without CO (2—1) lines are simply redshifted out of tunable 
range of the 1mm receivers and thus cannot be observed. We also calculated the 3a upper 
limits (about 10^*^ Kkms^^ pc^) for the two non-detections by assuming a line width (Full 
Width at Zero Intensity, hereafter FWZI) of 600 km s~^. The reasons for non-detections (for 
IRAS F18216+6419 and IRAS F02054-F0835) may be that the integration time is insufficient 
because these two objects have the highest redshifts in our sample. 

The CO (1-0) and CO (2-1) line widths (FWZI), integrated line luminosity, FIR 
luminosity and LFiR/L|-,Q(--,^_g^ for all our sample objects, including 6 IR QSOs from the 
literature, are listed in Table 2. From this table, we can see that the CO luminosities {L'qq) 
for the 23 CO detected IR QSO hosts (17 of which are from our sample) are several times 
10^ to a few times 10^*^ K km s^^ pc^, comparable to those for ULIRGs with the largest 
molecular gas reservoirs in the local universe. The upper limits of the two non-detections 
are also consistent with this mass range. 

The CO (2— 1)/C0 (1—0) line luminosity ratio (^21) ranges from 0.4 to 1.2 with mean 
value of ~ 0.8 in nine CO (2—1) detected IR QSOs. Four IR QSOs out of nine are found 
to have > 1 with uncertainty of ~ 15%, suggesting that the molecular gas is essentially 



- 8 - 



thermalized and has a point-like distribution compared to the 12" CO (2—1) beam (corre- 
sponding to 24kpc for the median redshift of CO (2—1) detected IR QSOs, z ~ 0.11). Three 
IR QSOs have r^x ~ 0.7 with uncertainty of ~ 12%, a typical average value for ULIRGs 
(Radford, Downes & Solomon 1991), which implies they, similar to ULIRGs, have modest 
sub-thermal molecular gas properties. Furthermore, if they are indeed sub-thermally excited 
with such line ratios, then these sources again have point-like molecular gas distribution com- 
pared to the 12" beam. Only 2 IR QSOs have really low r^x ratios (~ 0.5) with uncertainty 
of ~ 10%, which may be due to their sub-thermalized molecular gas properties and/or their 
extended spatial structure (Garay, Mardones & Mirabel 1993) of quite a few arc-seconds. 
From our limited CO (2—1) data, we infer that the molecular gas distributions in most (7 
out of 9) IR QSOs are probably concentrated in compact nuclear regions. However, further 
high spatial resolution observations, e.g., Aravena et al. (2011), are needed to reveal the 
spatial distribution and kinematics of the molecular gas in IR QSOs. 



4.2. Comparison with Local ULIRGs 

Fig. [3]compares the CO properties of IR QSO hosts (upper panels) with those of ULIRGs 
(lower panels). Specifically, the histograms of the CO luminosity ^qq{x-q)) (^^^^ panels), 
the FIR to CO luminosity ratio Lfir/I^co (commonly used as an indicator of the star for- 
mation efficiency [SFE], middle panels) and CO (1—0) line width (FWHM, right panels) are 
compared. Median and mean values are labelled in the upper left corner of the panels. 

From Fig. [3] we can see that the median value of for IR QSO hosts is similar 

to that of ULIRGs. Hence, the molecular gas in IR QSO hosts is as abundant as in ULIRGs 
which can easily provide the fuel needed for their ongoing starbursts. In addition, the median 
values of Lfir/I^co (SFE) and CO line width for IR QSO hosts are also similar to those of 
ULIRGs at the 3cr level. IR QSO hosts appear to have similar SFRs, velocity spread in CO, 
and SFEs as ULIRGs, likely though ULIRGs have evolved into the QSO episode manifested 
as IR QSOs. 

It is widely accepted that the FIR emission of ULIRGs is mainly from star formation (e.g. 
Yun et al. 2001; Gao & Solomon 2004b) and most molecular gas in ULIRGs is distributed 
in centrally rotating nuclear disks or rings as revealed by CO images (Downes & Solomon 
1998). The similarities of IR QSO hosts and ULIRGs in their CO properties suggest that 
their FIR luminosities are both dominated by dust heated by starbursts and even their cold 
molecular gas may have similar spatial distributions. In fact, the highest resolution CO 
interferometer observations for two IR QSOs, Mrk 231 and I Zw I show nuclear disks and 
rings. For Mrk 231, there is an inner disk of diameter 1.2" (radius ~ 0.5 kpc) and an outer 
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disk with diameter of 3" (radius ~ 1.25 kpc, Downes & Solomon 1998). For I Zw I, there is 
a ring-hke structure of radius 1.2 kpc (Staguhn et al. 2004). Our available line ratios of 
are also consistent with the CO emission in IR QSOs being mostly concentrated in nuclear 
regions. If most IR QSO hosts indeed have similar CO properties to those of ULIRGs, it 
follows that the CO-to-H2 conversion factor from L'qq to cold molecular gas mass for IR 
QSOs should be the same for ULIRGs, i.e. aco = O.SMq (K kms~^ pc^)~^. Therefore, the 
cold molecular gas properties, such as cold molecular gas mass, SFE, and possibly the spatial 
distribution in IR QSOs, are indistinguishable from those of ULIRGs. 

Taking into account the properties of IR QSOs in different wavebands, such as the 
optical morphology, extreme Fe II emissions in their optical spectra (Zheng et al. 2002), FIR 
excess emissions compared with classical QSOs (EOS) and similar 6.2/im PAH and [Nell] 
12.81/im emission strengths in IR QSO and ULIRGs (Cao et al. 2008), there is little doubt 
that there are ongoing massive starbursts in IR QSO hosts, as in ULIRGs. The timescale of 
starbursts in IR QSOs is about a few times 10^ yr, as inferred from Mgas/SFR by assuming 
that the FIR emission is predominately from star formation (see Table 3). 



4.3. Comparison with Local and High-redshift CO-detected QSOs 

Fig. m presents the Lpm versus -^co(i-o) i"6lation for almost all local and high-redshift 
QSOs with CO detections up to 2010. Throughout this paper, the red filled squares rep- 
resent the local IR QSOs, black open pentagrams denote the local PC and HE QSOs and 
blue crosses represent the high-redshift QSOs. It is clear from Fig. |4] that there is a tight 
correlation between Lpm and L'(^q^^_q^, spanning four and three orders of magnitude in Lpm 
and -Z^co(i-o)' respectively. The Spearman Rank-order (S-R) correlation analysis gives a 
coefficient of 0.94 with significance of > 99.99%. Such a strong correlation between Lpm 
and -^co(i-o) ^'^^ detected QSOs seems to suggest that the bulk of FIR emissions 

from all QSOs are from the same origin (star formation), since similar correlations are well 
established for star-forming galaxies at both low and high redshifts (e.g. Solomon & Vanden 
Bout 2005; Riechers 2011a). The best-fit power-law slope between LpiR and -^^co(i-o) ^P" 
proximately 1.4 regardless of the sample QSOs included. This is different from star-forming 
galaxies where the best-fit slope varies from 1 to 1.7 depending on how many ULIRGs are 
included (Gao & Solomon 2004b). 

The nonlinear correlation between LpiR and -^^co(i-o) ^'-'^ local star-forming galaxies, lu- 
minous starbursts, ULIRGs and high-redshift CO detected objects has been widely discussed 
and is likely due to different dense molecular gas fractions in different classes of objects (Gao 
& Solomon 2004b; Daddi et al. 2010). Our super-linear power-law slope (1.4) is consistent 
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with their conclusion. Furthermore, Fig. [5] shows the luminosity ratio Lfir/L^qj^j^_qj as func- 
tions of FIR luminosities Lfir (left panel) and CO (1-0) luminosities -^^co(i-o) (^ight panel). 
It is clear from Fig. [5] that there are tight correlations between Lfir/L^q(.-|^_q^ and Lfir, as 
well as between Lfir/L'(^q(.-^_p^ and L'fjQf^^_Qy Notice that unlike Fig. H], the vertical axis does 
not depend on the distance. The S-R correlation analysis gives a coefficient of 0.83 and 0.59, 
respectively, and both with significance level larger than 99.99%. Gao Sz Solomon (2004b) 
showed that Lpi^/L'^q^-^_q-^ is proportional to the dense molecular gas fraction as measured 
by Lhcn/Lco(i-o)' "where HCN is the most abundant high dipole- moment molecules. The 
correlations shown in Fig. |5] indicate that the dense molecular gas fraction might increase 
with increasing FIR luminosities and CO (1-0) luminosities (-^co(i-o)) detected 
QSOs. 

We also collected polycyclic aromatic hydrocarbon (PAH) 6.2/im emission data from 
the IRS on Spitzer (Houck et al. 2004) for our IR QSO sample and other CO detected QSOs 
(Lutz et al. 2008; Cao et al. 2008). Fig. |6]shows the correlation between Lq ^^ixh and -^co(i-o) 
for all sample QSOs with 6.2/im measurements including some upper limits. Since the 6.2/im 
PAH emission is believed to originate from star formation, the correlation shown in Fig. E] 
further confirms the star formation origin of the FIR emission for all QSOs. Netzer et al. 
(2007) reached the same conclusion based on their analysis of the SED of QSOs and ULIRGs 
observed by Spitzer/IR^. Our results further strengthen the conclusion: even in the QSO 
episode, the large amount of cold molecular gas reservoir still serves as the fuel to sustain 
star formation seen in the FIR emission. 

We now return to discuss Fig. |H A closer examination of this figure reveals that the 
correlation between Lfir and -Z^co(i-o) QSOs shown in Fig. H] is tighter than that for all 
star-forming galaxies, such as local spiral galaxies, ULIRGs and submm galaxies (SMGs, 
e.g.. Fig. 8 in Bertram et al. 2007). It appears to suggest that the cold molecular gas 
properties (distribution) in QSO hosts is not as complicated as those in starbursts, ULIRGs 
and/or SMGs. 

This point is supported by different spatial distributions of molecular gas in these ob- 
jects. High resolution CO observations for local IR QSOs (e.g., Mrk 231 and I Zw I, as 
described in §4.2) and high-redshift (sub)mm detected QSOs (e.g. APM 08279-1-5255, Clover- 
leaf QSO, J1148-I-5251) show that the molecular gas in these QSO hosts are concentrated 
in central disks/rings with radius on the order of 1 kpc (Solomon et al. 2003; Staguhn et 
al. 2004; Riechers et al. 2009; Walter et al. 2009). Bradford et al. (2009) reached the same 
conclusion from analyses of many molecular line species from the single-dish observations 
of the Cloverleaf. In contrast, a significant fraction of ULIRGs and SMGs are interacting 
or merging galaxy pairs with wide separations. In this case, a lot of the cold molecular 
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gas has not yet settled into nuclear rotating disks/rings, instead, the gas is still very ex- 
tended spatially and likely has different physical parameters from the molecular gas in the 
disks (Dinh-V-Trung et al. 2001; Greve et al. 2005; Tacconi et al. 2008; Riechers et al. 
2011c), which may have contributed to the larger scatters in the correlation between Lpm 
and L'(.Q(^_Q). 

Further evidence for different sizes is provided by recent high resolution observations 
of SMGs and lensed high-z QSOs (Casey et al. 2011; Riechers et al. 2011c, 2011d), which 
reveal extended cold molecular gas components in SMGs, but not in lensed QSOs. The 
tighter correlation between Lpm and -Z^co(i-o) ^'^^ QSOs, as shown in Fig. HI is consistent 
with different spatial distributions of molecular gas between QSOs and ULIRGs/SMGs. 

From Fig. HI one can also see that while most CO detected QSOs at high-redshift are 
more luminous, nevertheless the local and high-redshift CO detected QSOs form a contin- 
uous sequence. In particular, four high-redshift CO detected QSOs have CO luminosities 
comparable to local IR QSOs. The comparable CO and FIR luminosities imply similarities 
between these high-redshift QSOs and local IR QSOs in both cold molecular gas contents 
and SFRs. Due to sensitivity limits of current facilities, observations of FIR and CO faint 
QSOs are usually for gravitationally lensed systems: three of these four relatively faint 
high-redshift QSOs are lensed (SMM 04135, IRAS F10214 and Cloverleaf). Most recently, 
Riechers (2011a) reported CO detections for four additional lensed high-redshift QSOs with 
relatively low Lfir and -^co(i-o) Combined Array for Research in Millimeter-wave 

Astronomy (CARMA). Again, the intrinsic FIR and CO luminosities of these high-redshift 
QSOs are comparable to our IR QSOs. The similarities can be further verified when more 
observations of fainter (unlensed) high-redshift QSOs become available with the Atacama 
Large (sub-) Millimeter Array (ALMA) in the near future. 



4.4. the Evolution from Massive Starbursts to Luminous QSOs 

It is widely accepted that the growth of central SMBHs and their host spheroids are 
closely related due to the local SMBH-spheroid mass relation (Magorrain et al. 1998; Fer- 
rarese & Merritt 2000; Tremaine et al. 2002). This close connection is also mirrored in 
similar shapes in the cosmic history of black hole growth and star formation (e.g. Merloni 
et al. 2004). It is suggestive of a possible coeval evolution of central SMBHs and their 
spheroids. However, it is still unclear how or even whether spheroids of galaxies co-evolve 
with their central SMBHs at all stages. In the following, we examine more closely the rela- 
tion between the accretion rate of central SMBHs and SFRs in the host galaxies for all CO 
detected QSOs. 
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Fig. [7] shows the AGN-associated bolometric luminosities (Lagn) vs. -Z^co(i-o) 
CO detected QSOs. Note that the AGN-associated bolometric luminosities used in this 
paper are the integrated luminosities over the SEDs for the AGN components, as described 
in §3. From Fig. [TJ we can see clearly that as -^co(i-o) increases, Lagn also increases. 
The Spearman Rank-order correlation analysis gives a coefficient of 0.83 with significance 
> 99.99%. A least-square bisector fit for low-redshift CO detected PC and HE QSOs in 
Fig. [7] yields a best-fit slope of 1.4 ± 0.2. The correlation between Lagn and -^co(i-o) 
suggestive of some link between the cold molecular gas on ~ kpc scale disks/rings and the 
central SMBH accretion process on much smaller scales - it is likely that both nuclear star 
formation and SMBH accretion rely on the same molecular gas reservoir (Bonfield et al. 
2011; Diamond-Stanic & Rieke 2011). 

It is interesting that all IR QSOs and the few relatively faint high-redshift QSOs with 
Lagn less than ~ IO^^Lq are located on or below the best-ffi line, while most high-redshift 
QSOs with Lagn larger than ~ IO^^Lq are located above. Such a systematic difference in 
the locations on Fig. [7] for IR QSOs and most high-redshift optically bright QSOs could 
reflect some intrinsic difference in the properties between these two populations, although 
some selection effects may have also played a role. 

The Lagn/-^co ratio relates to the accretion efficiency to central SMBH, as indicated 
by Macc/Mgas- Fig. [8] compares the histograms of Lagn/L'^q for CO detected local IR QSOs, 
PG-I-HE QSOs and high-redshift QSOs (top, middle and bottom panels respectively). A 
comparison of the top and middle panels shows that the median value of Lagn/-^co for local 
PG+HE QSOs is larger than that for IR QSOs. Since a large fraction of optically bright 
PC QSOs are CO non-detected (or not observed), the real median value of Lagn/Lqq for 
classical QSOs should be even higher than the plotted one. The difference of the Lagn/-^co 
ratios between IR QSOs and the brightest high-redshift QSOs is even more dramatic, about 
one order of magnitude. 

In order to understand why the Lagk/Lqq ratios for local PG+HE QSOs and high- 
redshift QSOs are much higher than those of IR QSOs, we plot the histograms of LfirZ-^^agn 
in Fig. inifor these three classes of objects. For local classical QSOs and the brightest high- 
redshift QSOs, the FIR luminosities are ^10% of their central AGN bolometric luminosities, 
so the main energy output for these QSOs is dominated by central AGNs, rather than 
starbursts. In contrast, the median Lpm/ijAGN value for IR QSOs is ~0.6, indicating that 
the FIR luminosities of most IR QSOs are comparable to their central AGN bolometric 
luminosities. 

The ratio of central SMBH accretion rate to the SFR in their host galaxies (Macc/SFR) 
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can be estimated by Lagn/-Z^fir based on the following equation: 

Ma,c/SFR = 2.1Lagn/^fir x 10-^ (2) 

assuming an accretion efficiency rj = 0.1 (see Hao et al. 2008). Although the uncertainties 
are large due to a number of assumptions on determining these key parameters, we can 
nevertheless roughly estimate whether SMBHs and their host galaxies co-evolve based on 
the Lfir/Lagn ratio. 

As can be seen from Fig.[9l the median values of Lfir/-^vagn are about 0.58, 0.16 and 0.13 
for IR QSOs, local PG+HE QSOs and high-redshift QSOs, respectively. The corresponding 
median values of Macc/SFR are 3.6 xlO-^ for IR QSOs, and 1-2 xlO-^ for local PG+HE 
QSOs and high-redshift brightest QSOs (see also Coppin 2009). So the Macc/SFR value for 
IR QSOs is comparable to that of local MBu/Msph value, ^ 1.4 x 10"^ (Haring & Rix 2004), 
implying their central SMBHs and host spheroids grow in line with the local relation. In 
contrast, the Macc/SFR values for local classical QSOs and high-redshift bright QSOs are 
almost an order of magnitude larger than the local M-QYi/Msp^ value. For these objects, the 
central SMBH masses increase much faster than their host spheroids compared to the local 
scaling. It is interesting to note that the four relatively faint high-redshift CO detected QSOs 
(SMM 04135, F 10214, Cloverleaf and RXJ 124913) have high Lfir/Lagn values that are 
close to or even larger than those of local IR QSOs, thus the growth of their black hole may 
be slower than that of high-luminosity QSOs. 

According to the scenario ffist pointed out by Sanders et al. (1988a, b), IR QSOs will 
eventually evolve toward classical QSOs. The different efficiencies of gas accretion to central 
SMBH (Macc/Mgas) between the IR QSOs and local classical QSOs and the high-redshift 
brightest QSOs then ffi well into a picture where the QSO luminosity reaches the maximum 
when the star formation at the central ~ kpc decreases to a few tenths of its peak value 
and no longer dominates the energy output (as indicated by the Lfir/I^agn value). Such an 
evolutionary path has also been discussed for AGNs in SMGs by Alexander (2009) and is 
broadly consistent with the simulation results of Hopkins (2011). 

We can also obtain more quantitative estimates of the star formation timescales and 
black hole masses. Table 3 lists the median values of SFR, Mace, Mbh, cold molecular 
gas mass Mgas and dynamical mass Mdyn for IR QSOs, PG+HE QSOs and high-redshift 
QSOs. The SFRs and accretion rates are estimated from Lfir and Lagn following Hao et al. 
(2008), the SMBH masses are taken from the literature, while gas masses are estimated from 
-^co(i-o) dynamical masses are estimated from the stellar CO absorption line widths in 
the H-band (Dasyra et al. 2006) for local QSOs and CO line widths for high-redshift QSOs. 
Such dynamical masses should be treated as upper limits of the stellar mass of QSO's host 
galaxies (e.g. Coppin et al. 2008). 
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First, we notice that the timescale for the gas consumption is on the scale of a few- 
times 10^ yr for all types of objects. Second, the MBn/A^dyn for PG QSOs is 10~^, similar to 
the value for local galaxies, again indicating that PG QSOs follow the local scaling relation 
(Dasyra et al. 2006). However, the MBn/^dyn values arc quite different, 4.5 x 10~^ and 10^^ 
for IR QSOs and the brightest high-redshift CO detected QSOs, respectively. The SMBHs 
in IR QSOs appear to be below the local Mbh and Mgph relation, whereas the SMBHs in the 
brightest high-redshift CO detected QSOs are much larger than that predicted by the local 
relation. It suggests that while star formation and AGN activities are intimately connected, 
they are not necessarily synchronized, as also suggested by numerical simulations (Hopkins 
2011). In particular, IR QSOs appear to be in a transition stage with a timescale of a few 
times 10^ yr; the rapid accretion (with Mace ~ yr~^) can grow the mass of a black hole 
up to a few times IO^Mq after this phase. Henceforth it may follow the Mbh vs. Mgph 
relation seen in local galaxies. 



5. SUMMARY 

In this paper we report the detections of cold molecular gas in 17 out of 19 infrared 
ultraluminous QSO (IR QSOs) hosts observed using the IRAM 30m telescope. Including 
six additional IR QSOs with CO detections from the literature, our IR QSO sample consists 
of 25 objects with 23 detections. The gas reservoir in these objects is found to be several 
times IO^Mq to a few times lO^^M©. Comparisons with local ULIRGs indicate that the 
cold molecular gas properties, such as the cold molecular gas mass content, star formation 
efficiency, and the CO (1-0) line widths of IR QSO hosts are similar to those of ULIRGs. 
These results suggest that sufficient amount of cold molecular gas exists to sustain massive 
starbursts even in the ultraluminous IR QSOs phase. 

We also compared the IR QSO properties with several other QSO samples with molec- 
ular gas detections at both low and high redshifts. Our main findings are summarized as 
follows. 

1. There exists a tight correlation between Lpm and -C'co(i-o) QSOs, spanning four 

and three orders of magnitude in LpiR and L'^q(^^_q^, respectively. The 6.2 fim PAH 
luminosity (L6.2^m) and L'^^q^^ ^^ are also closely correlated for all sample QSOs. It 
seems to confirm the results based on Spitzer/IRS observations that the FIR emissions 
of all QSOs are mainly from star formation process rather than AGN. From this cor- 
relation as well as our limited CO line ratio measurements for IR QSOs and available 
CO images of QSOs, we speculate that the geometry of cold molecular gas for most 
QSOs at both high-redshift and low-redshift could be in a disk/ring on ~ kpc scale. 
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similar to those seen in some local ULIRGs. 

2. The AGN-associated bolometric luminosities of all QSOs increase as Lqq increases, 
implying a possible link between the cold molecular gas on ~ kpc scale (disk/ring) and 
the central black hole accretion process. Thus it is likely that both star formation and 
central black hole accretion draw from the same cold molecular gas reservoir. 

3. The Macc/SFR values for IR QSOs and a few high-redshift, relatively faint QSOs are 
comparable to the local MsH/Msph value. These QSOs might be in the transition 
stage from gas-rich galaxy mergers to QSOs then to elliptical galaxies, exhibiting both 
high SFR and high accretion rates. However, the local MeH/Msph relation could not be 
established in this short transition phase. If the black hole continues to grow vigorously 
after this transition period, then the MBn/Msph relation may be established afterwards. 
On the other hand, for both local and very bright high-redshift QSOs, the black hole 
appears to grow much faster than the spheroids. It remains to be seen how the Mbh 
vs. Mgph evolves as a function of cosmic time. 

This work is based on observations carried out with the IRAM 30m telescope. IRAM is 
supported by INSU/CNRS (Prance), MPG (Germany) and IGN (Spain). We thank Drs. J. 
S. Huang and J. M. Wang for advice and helpful discussions. We also thank Drs. Thomas 
Bertram and Lutz Wisotzki for kindly providing us with the unpublished Bj magnitudes for 
HE QSOs, and Dr. Ran Wang for providing us data for high-redshift QSOs. This project 
is supported by the NSF of China 10833006, 10973011, 11003015, 11173059 and 973 project 
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Fig. 1. — True color images of 12 IR QSOs (out of a total of 25), constructed from SDSS 
g, r, and i images using color-preserving nonlinear stretches (Lupton et al. 2004). The box 
size of each image is 80 kpc. 
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Fig. 2. — Spectra of CO (1-0) (solid lines) and CO (2-1) (dotted lines, divided by a factor 
of 4 for 9 targets) emission for our IR QSO sample. The spectra are smoothed to 23 km 
channel width for display. Only marginal detections (3cr) in F12265 and F12134, and 
non-detection in F02054 and F18216; for these objects, the spectra are smoothed to 46 km 
s~^. The main-beam temperature can be converted to the flux density using S/Tnib=4.95 
Jy/K. 
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Fig. 3. — Histograms of CO (1-0) luminosity Lco(i-o)' ^^^^ formation efficiency (ratio of 
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IR QSO hosts (upper panels) and ULIRGs (lower panels). The median and mean values are 
indicated in the top-left corner of each panel. 
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Fig. 4. — LpiR vs. -Z^co(i-o) almost all CO detected IR QSOs (red filled squares), PG+HE 
QSOs (black open pentagrams) and high-redshift QSOs (blue crosses). The samples are 
described in §2. A least-square bisector best-fit line obtained for all objects is shown; the 
best-fit power-law slope is 1.4 ± 0.1. 
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Fig. 5. — The star formation efficiency Lfir/L[^q^-^_qj vs. Lfir (left panel) and Lfir/L[^q^^_q^ 
vs. -Z^co(i-o) ("ght panel) for three samples of QSOs (see §4.2). A Spearman Rank-order 
correlation analysis gives significances of > 99.99% for both correlations. Notice that the 
vertical axis is independent of distance. 
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Fig. 6. — PAH 6.2/im luminosity L6.2^m vs. -^co(i-o) 

for PAH and CO detected IR QSOs, 
PG+HE QSOs and high-redshift QSOs. The PAH luminosities of low-redsliift and liigh- 
redshift QSOs are taken from the Quasars and ULIRGs Evolution Study (see Schweitzer et 
al. 2006; Netzer et al. 2007) and Lutz et al. (2008) respectively, whereas PAH luminosities 
of IR QSOs are from Cao et al. (2008). Upper limits are plotted for some objects. A least- 
square bisector best-fit line is obtained for all objects, including those with upper limits. 
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Fig. 7. — AGN-associated bolometric luminosity (Lagn) vs. Lco(i-o) three samples of 
QSOs. A least-square bisector best-fit was performed for only low-redsliift CO detected 
PG+HE QSOs; the best-fit power-law slope is 1.4 ± 0.2. 
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Fig. 8. — Histogram of the Lagn to L'qq ratio for IR QSOs (top panel), PG+HE QSOs 
(middle panel) and high-redshift QSOs (bottom panel). The median and mean values are 
indicated in the top right corner of each panel. Notice that a large fraction of optically 
bright PG QSOs are not detected in CO, thus the true median value of -^^agn/-^co t^-king 
into account the non- detections will be larger than that shown in the middle panel. 
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Fig. 9. — Histogram of I^firZ-^agn for IR QSOs (top panel), PG+HE QSOs (middle panel) 
and high-redshift QSOs (bottom panel). The median and mean values are indicated in the 
top right corner of each panel. 



Table 1. CO observations of IR QSOs with the IRAM 30m telescope 



Source 


R.A. 


Decl. 


^opt 


^CO 




Transition 


*f„t 




(J2000.0) 


(J2000.0) 






(Mpc) 




(min) 



Detections 



F02066+470B 


02 


09 


45. 


,8 


+47 


19 


43. 


2 


0. 


132 


0. 


132 


620 


CO(l-J-O) 
CO(2->.l) 


47 
47 


F04415+1215 


04 


44 


28. 


8 


+ 12 


21 


13, 


1 


0, 


089 





090 


407 


CO(l^O) 
CO(2^1) 


85 
66 


IR()6269-0543 


06 


29 


24. 


7 


-05 


45 


26, 


,0 


0. 


117 





117 


545 


CO(l^O) 
CO(2^1) 


97 
78 


F07599+6508 


08 


04 


30. 


,4 


+64 


59 


53, 


,3 


0. 


148 


0, 


,149 


703 


CO(l^O) 
CO(2-!.l) 


304 
219 


F10026+4347 


10 


06 


41. 


,8 


+43 


32 


41. 


,6 


0. 


178 


0. 


,179 


861 


CO(l->-0) 


99 


F11119+3257 


11 


14 


38. 


,9 


+32 


41 


33. 


,0 


0. 


189 


0. 


,190 


920 


CO(l^O) 


94 


Z11B98-0112 


12 


02 


26. 


,6 


-01 


29 


15. 


,3 


0. 


161 


0. 


,151 


718 


CO(l^O) 
CO(2^1) 


115 

72 


F12265 + 0219 


12 


29 


06. 


6 


+02 


03 


09, 


,0 


0. 


158 





159 


755 


CO(l^O) 


150 


F12134+5459 


12 


15 


49. 


,3 


+54 


42 


24, 


,6 


0. 


150 


0, 


151 


713 


CO(l^O) 


140 


F13218+0552 


13 


24 


19. 


9 


+06 


37 


05. 


,0 


0. 


,206 


0. 


,204 


1007 


CO(I-J-O) 


113 


F13342+3932 


13 


36 


24. 


,0 


+39 


17 


32. 


,2 


0. 


179 


0. 


,180 


866 


CO(l->-0) 


190 


F16069+1808 


16 


09 


13. 


,7 


+ 17 


57 


11. 


,0 


0. 


171 


0. 


,171 


824 


CO(1->0) 


115 


F16462-0460 


16 


48 


66. 


,8 


-04 


69 


33. 


,5 


0. 


101 


0. 


,100 


465 


CO{1^0) 
CO{2^1) 


143 
91 


F20036-1547 


20 


06 


31. 


9 


-15 


39 


05 


8 


0. 


193 





192 


942 


CO(l^O) 


98 


F21219-1757 


21 


24 


41. 


,6 


-17 


44 


45, 


,3 


0. 


113 


0. 


113 


525 


CO(l^O) 
CO(2-).l) 


198 
137 


F22464-1744 


22 


48 


04. 


1 


-17 


28 


28. 


,5 


0. 


117 


0. 


,118 


545 


CO(l->-0) 

CO(2-+l) 


66 

66 


F23411+0228 


23 


43 


39. 


,7 


+02 


45 


05. 


,7 


0. 


,091 


0. 


,092 


416 


CO(l^O) 
CO(2^1) 


113 

94 



Non-detections 



F02054+0835 02 08 06.8 +08 50 05.2 0.345 ■■■ 1826 CO(l-)-O) 38 

F18216+6419 18 21 57.3 +64 20 36.0 0.297 ■ ■ ■ 1535 CO(l-)-O) 58 



Note. — Units of right ascension and declination are hours, minutes, seconds, and degrees, 
arcminutes, arcseconds respectively. 

^CO redshift was determined by computing its flux-weighted redshift, zqq = ^I{z)z/^I{z) 
(see, e.g., Greve et al. 2005). 

^ Dj^ is the luminosity distance. 

'^Totai usable on-source integration time. 



Table 2. CO emission line and infrared properties of the IR QSO sample 



Source 


Transition 




^^FWZI 








logiFIR 


ipia/ico 






(km 


(km s-1) 


(K km s"'^) 


(Jy km s-1) 


(lO^K km s-^ pc^) 


(logi0) 




Detections 


F02065+470B 


GO(l^O) 


145±12 


300 


1.84+0.17 


9.11+0.84 


7.60 


11.886 


101 




GO(2^1) 


172±10 


330 


6.19+0.29 


26.69+1.44 


6.36 






F04415+1216 


GO(l->-0) 


112±8 


230 


1.49+0.12 


7.38+0.59 


2.76 


12.011 


373 




GO(2->.l) 


137±16 


280 


2.48+0.29 


12.28+1.44 


1.14 






F06269-0543 


CO(l^O) 


258+25 


465 


2.28+0.22 


11.29+1.09 


7.36 


12.103 


172 




CO(2^1) 


320+22 


540 


11.00+0.68 


54.45+3.37 


8.87 






F07599+6508 


CO(l^()) 


376+14 


770 


2.36+0.12 


11.68+0.59 


12.33 


12.074 


96 




CO{2^1) 


486±44 


840 


7.36+0.78 


36.43+3.86 


9.60 






F10026+4347 


CO{1^0) 


343±92 


680 


0.81+0.16 


4.01+0.79 


6.20 


11.830 


109 


F11119+3257 


CO(l^O) 


285+36 


470 


1.31+0.16 


6.48+0.74 


11.33 


12.273 


165 


Z11598-0112 


CO(l-)-O) 


207+16 


680 


2.43+0.18 


12.03+0.89 


13.24 


12.260 


137 




CO{2^1) 


234+18 


580 


9.71+0.70 


48.06+3.46 


13.19 






F12265+0219-'' 


CO{1^0) 


313+86 


480 


4.88+1.41 


24.16+6.98 


29.11 


12.266 


63 


F12134+5459^ 


CO{1^0) 


359+84 


525 


0.56+0.16 


2.77+0.79 


2.94 


11.792 


211 


F13218+()r)r)2 


CO(l^()) 


299+64 


430 


0.93 + 0.18 


4.60+0.89 


9.51 


12.175 


157 


F13342+3932 


CO{1->0) 


333+49 


520 


1.99+0.11 


9.85+0.54 


15.44 


12.125 


86 


F15069+1808 


CO(l->-0) 


287±B0 


440 


1.96+0.14 


9.65+0.69 


13.71 


11.876 


66 


F15462-0450 


CO(l->-0) 


208+18 


490 


1.41+0.13 


6.98+0.64 


3.36 


11.964 


267 




CO(2->l) 


177+18 


370 


4.25+0.32 


21.04+1.58 


2.64 






F20036-1547 


CO{1^0) 


263+23 


420 


1.91+0.18 


9.45+0.89 


17.20 


12.338 


127 


F21219-1757 


CO{1^0) 


275+82 


400 


1.11+0.09 


5.49+0.45 


3.34 


11.631 


128 




CO(2^1) 


243+28 


380 


4.19+0.19 


20.74+0.94 


3.14 






F224B4-1744 


CO(l-)-O) 


186+13 


380 


2.27+0.19 


11.23+0.94 


7.31 


11.648 


48 




CO(2->.l) 


201+7 


460 


9.30+0.36 


46.04+1.73 


7.49 






F23411+0228 


CO(1->0) 


218+22 


430 


1.44+0.16 


7.13+0.79 


2.78 


11.742 


199 




CO(2->.l) 


229+21 


390 


2.98+0.30 


14.76+1.48 


1.43 






Non-detections 


F()2()r)4 + ()835 


CO(l^O) 






< 0.44 


< 2.17 


< 12.91 


12.655 


> 350 


F18216+6419f 


CO{1->0) 






< 0.70 


< 3.46 


< 15.08 


12.694 


> 328 


Literature sources 


PG 0060+124'' 


GO(l-!-0) 


370 




6.0+0.2 


30+1.1 


6.12 


11.417 


61 


Mrk 1014'' 


GO(l->-0) 


270 




1.1+0.1 


6.6+0.6 


6.94 


12.279 


274 


Mrk 231' 


CO(l^O) 


230 




22.0 


99.0 


8.98 


12.185 


171 


PG 1613+658'' 


CO{l^()) 


400 




1.6+0.1 


8.0+0.6 


6.23 


11.595 


63 


3C48J 


CO(l^O) 


330 






1.9+0.2 


18.49 


12.666 


251 


PG1700+518'' 


GO(l-!-0) 


260 




0.79+0.14 


3.9+0.7 


16.07 


12.186 


96 



*The full width at half maximum (FWHM) obtained from Gaussian fit. 

^Thc full width at zero intensity (FWZI) of the CO emission line width. 

^The total CO line intensity /co obtained by integrating over the full velocity range, the errors were calculated using equation (1) in 

Matthews & Gao (2001). 

^■S'cO^" the CO flux obtained using the conversion factor S/T^^— 4.95 Ji^ ~ ^ 

°L^O = 2.45 X 10^ ( jjf kmt-1 ) {v^y iJ^^il + 2)"^)IK km s'^ pc^] (see, e.g., Solomon et al. 1997). 
^Marginally detected with a signal-to-noise ratio of ~3. 

SAravena et al. (2011) report the GO (1-0) detection for PG1821+643 by CARMA. The detected CO(l-O) luminosity is about one half of 
the 3cT upper limit shown here. 

''Evans et al. (2006) 

'Solomon et al. (1997); Kim & Sanders (1998) 
■iKrips et al. (2006); Heckman et al. (1992) 
■•Evans et al. (2009) 



Table 3. Median values of physical parameters for IR QSOs, PG+HE QSOs and 

high-redshift QSOs 



Name 


^FIR 

(10^^ is) 


SFR 
(Moyr-i) 


Lagn 


M 

(Moyr-i) 


(10^ Mo) 


gas 
(109 Mq) 


^dyn 
(10" Mq) 


IR QSO 


1.2 


391 


1.4 


0.9 


5 


6 


0.9 


PG+HE QSO 


0.05 


16 


0.2 


0.13 


20 


0.8 


2 


high-redshift QSO 


7 


2300 


100 


70 


200 


20 


2.1 



Note. — The median values of parameters shown in this table are not obtained for all our sample objects 
due to incomplete information for some. 

'^The black hole masses for local IR QSOs and PG QSOs are from Hao ct al. (2005); for high-redshift 
QSOs, the median value of black hole masses are taken from Coppin et al. (2008) for QSOs at z ~ 2. Prom 
the most recent measurement by De Rosa et al. (2011), the black hole masses for QSOs a,t A < z < 6.5 are 

^10^ Mq, smaller than those at lower redshift. Notice that the median value of the black hole masses for the 
local classical QSOs are only based on PG QSOs due to the lack of Mbh information for HE QSOs. 

''The molecular gas mass for all QSOs is estimated from -f'co(i-o) '^ith a conversion factor of aco = 
O.8M0(K kms-i pc^)-^ 

'^The dynamical masses of IR QSOs and PG QSOs are from Dasyra et al. (2006, 2007). Notice that we 
assume the median value of dynamical mass for IR QSOs is the same as that of local ULIRGs. The dynamical 
masses for high-redshift QSOs are estimated by the median value of their CO line width (see Coppin et al. 
2008). 



